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Abstract. Alternative gravitations of Milgrom (MOND), Moffat (MOG), and CDM scenarios all simulate rotation 
curves of spirals with reasonable details. They display significant disparities however in predicting the stellar mass- 
to- light (Mt/L) ratios of the galaxies. We maintain this feature could serve as a distinguishing factor between 
different alternative theories. We analyze the rotation curves of 46 low- and high-surface brightness galaxies and 
compare the resulting Af, /Ls with the predictions of the Stellar Population Synthesis (SPS) scheme. The color- 
Mt/L correlation obtained for MOND is consistent with predictions of SPS models. MOG does not show this 
consistency, and the K'h/Ls of CDM model shows large dispersions. Furthermore, M-t/L ratios of MOND with 
Bekenstein interpolating function favor Kroupa's initial mass function (IMF) of the SPS scheme, while those of 
MOND with standard and simple interpolating functions are consistent with Salpeter's IMF. Here is another 
indication to differentiate between different IMFs that are used in SPS context. 
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1. Introduction 

The gravitational force of the observable mass of large 
astronomical systems, galaxies, clusters of galaxies, or 
for that matter, the universe in general, is not suffi- 
ciently strong to explain the observed dynamics of the 
systems. To resolve the dilemma, one main school of 
investigators has resorted to dark matter/dark energy 
scenarios. In spite of extensive efforts, however, no one 
has so far reported a direct identification of the hy- 
pothesized dark entity through non-gravitational inter- 
actions with the observable matter. This lack of di- 
rect identification has inspired an equally intensive ef- 
fort to contemplate alternative theories of gravitation. 
The Modified Newtonian Dynamics (MOND) of Milgrom 
• (1983) and of Bekenstein (2004), the Modified Gravity 
(MOG) of Moffat (2005), the Nonlocal Nonlinear gravity 
of Sobouti (2008a, b, 2009), and varieties of f{R) grav- 
ities dCapoz ziello 2002*, Capozziello e t al. 2006 k. 2007[ 
ICarroll et aT 2004nSobouti 2007. .Sobouti et al. 2009P fall 



in this category. 

Rotation curves of spiral galaxies as measured by the 
21 cm line of HI often extend well beyond the optical disks 
of the galaxies and provide a valuable body of data to de- 
termine the radial dependency of the gravitational forces 
in galactic scales. In this paper we construct rotation 
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curves of a large sample of galaxies from the distribution 
of their detectable matter through three different gravity 
models, MOND, MOG, and Newtonian gravity plus cold 
dark matter (CDM) halos. At first glance all three mod- 
els seem to reproduce the observed data with reasonable 
detail. On a deeper examination, however, we find signifi- 
cant disparities in their predictions of stellar mass-to-light, 
M*/L, ratios. To differentiate between the models we re- 
sort to stellar population synthesis, SPS, analysis and the 
color- M,/i correlation predicted therein through various 
initial mass functions (IMF). There is the possibility to 
use this feature to discriminate between different gravity 
models and different IMFs. 

The paper is organized as follows: In Sect. 2 we give 
a brief review of the different gravity models used in our 
analysis. In Sect. 3 we describe our galaxy sample. Fits 
to the observed rotation curves are discussed in Sect. 4. 
Numerical results and brief concluding remarks are given 
in Sects. 5 and 6. 

2. Alternative gravity models 

In this section we review the basic tenets of two alternative 
gravities as well as the Newtonian gravity plus CDM halos. 
We present the end formulas that we will use in the study 
of the dynamics of galaxies. All three accommodate the 
two main asymptotic features of the rotation curves of 
spirals: the slow non-Keplerian decline of the curves at 
large distances from the galaxy, and the Tully-Fisher (TF) 
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relation, which is the approximate proportionality of the 
asymptotic speed of an orbiting object to the fourth root 
of the mass of the galaxy ( |Tully fc Fisher IQ??] ). 

Actually there is much debate on how fast or slow 
the rotation curves decline, if at all (|Persic et al. 19961 
ISalucci et al. 20071 IGentile 2008]) . There are also refine- 
ments and redefinitions to the Tully-Fisher relation. 
McGaugh (2005) prefers to use the total baryonic (stellar 
-|- gaseous) mass in TP relation to accommodate the gas- 
rich galaxies. See also Stark et al. (2009), for a calibration 
of the baryonic TF relation with the help of gas dominated 
galaxies. Nonetheless, both assumptions are adequate ap- 
proximations to the observed data and will be employed 
in this paper. 

2.1. Modified Newtonian Dynamics of Milgrom, 
MONO 

Based on observations of galactic rotation curves, Milgrom 
(1983) argues that the Newtonian dynamics is not vi- 
able below a certain universal acceleration, oq ~ 1.2 x 
10~^°m/sec^. To comply with the Tully-Fisher relation he 
modifies the law of motion to have an asymptotic acceler- 
ation proportional to the square root of the Newtonian 
acceleration. The MOND acceleration, gmond, and the 
Newtonian one, qn, are connected through Eq. (P) below 



Qmond I 9mond 
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ao 



9N 



(1) 



where fi{x) is an interpolating function for transition from 
the Newtonian to the MONDian regime. It runs smoothly 
from fj.{x) = X for cc << 1 to ^{x) = 1 for a; >> 1. Here 
we adopt three functions commonly used, the standard 
interpolating function of Bekenstein & Milgrom (1984): 

/^i(-) = (TTW^' 

the simpler function of Famaey & Binney (2005), 



and Bekenstein's interpolating function ([Bekenstein 2004t 

(l + 4x)i/2_i 



^3(2^) 



(1 + 4a;)i/2 + r 



(4) 



Hereafter, the analysis using /ii, /i2 and /^s will be refereed 
to as MONDl, M0ND2 and M0ND3, respectively. 

2.2. Modified Gravity of Moffat, MOG 

Modified Gravity of Moffat consists of three theories 
of gravity: the nonsymmetric gravity theory (NGT), 
the metric-skew-tensor gravity (MSTG) theory, and the 
scalar-tensor- vector gravity (STVG). They rely on the ex- 
istence of a massive vector field universally coupled to 
matter. Moffat maintains that MOG explains the rota- 
tion curves of galaxies, clusters of galaxies, and cosmolog- 
ical issues without resorting to dark matter (|Moffat 19951 



IMoffat 20051 IMotfat 20061 IMotfat fc Toth 20091) . Good 
fits to astrophysical and cosmological data have been ob- 
tained with his recent version of STVG. One notable fea- 
ture of NGT, MSTG, and STVG is that the modified ac- 
celeration at weak gravitational fields has a Yukawa-type 
addition to the Newtonian acceleration. In the weak field 
limit, STVG, NGT, and MSTG produce similar resuhs. 
The recipe for the gravitational force of a spherically dis- 
tributed mass, M(r), is (|Moffat 2006P 
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where Gjv is the Newtonian gravitational constant, M{r) 
is the baryonic mass inside the radius r, and a(r) = 
[Mo/M(r)]i/2^ The parameters Mq 

or To determine the 
coupling strength of the vector field to the baryonic mat- 
ter and to the range of the force, respectively. They are 
not universal constants and vary with the size of the sys- 
tems (jBrownstein fc Mofl: at ^0061 [Haghi fc Rahvar 2010"] ). 
In galactic scales, they are determined by analyzing the 
best fit of the theory to the rotation curves of LSB and 
HSB galaxies. For normal size galaxies, they are reported 
as Mo = 9.6 X IO^A/q and tq = 13.9 kpc, and for 
dwarf galaxies, as Mq = 2.4 x IO^^A/q and rg = 9.7 
kpc ( Moffat 2QQ6j . An empirical fitting of Mq versus tq 
for a wide range of spherically symmetric systems, from 
solar size to clusters of galaxies is depicted in Fig. 2 of 
Brownstein & Moffat (2006). The MOG gravitation tends 
to the Newtonian one as Afo and tq — !■ 00. 

2.3. Newtonian Gravity plus Cold Dark Matter, COM 

In this scenario, gravitation is Newtonian. To account for 
the nonclassical behavior of the rotation curves one adds 
a spherically symmetric dark halo to the galaxy. Here, we 
consider a NFW halo with the density distribution 
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and the gravitational acceleration 



gNFW = 47rGpsrs( — ) 
r 



In 1 



(l + r/r,) 



(6) 



where and ps are the characteristic radius and density 
of the distribution ( [Navarro, Frenk fc White 1996^ . The 
NFW density comes from numerical simulations of ACDM 
theory in the framework of structure formation. There, 
one also finds that these parameters are correlated to each 
other as in Eqs. (I7|and|8|) below, leaving only one free pa- 
rameter to characterize the halo (see IBuUock et al. 200T1 
IWechsler et al. 20021 and INeto et al. 20071 for details). 
Thus, 



A 



3 ln{l + c) -c/(l + c) 



(7) 



Hasani et al.: SPS scheme, a discriminant between gravity models 



3 



c= 13.6 



-0.13 



0.46 



kpc, (8) 



where pc is the critical density of the Universe and 
A = 200 is the virial overdensity at redshift z=0 
dBryan fc Norman 19981 ) . 



3. Observational data 

There are diverse morphological types of galaxies with 
diverse shapes and sizes to their rotation curves. Our 
sample, a collection of 46 galaxies taken from Sanders 
(1996), McGaugh & de Blok (1998), Sanders & Verheijen 
(1998), and Begeman (1991), accommodates these diversi- 
ties. Members of the sample have well measured rotational 
speeds and accurate surface photometry. They are listed 
in Table [Hand shown in Figs. [1]-E1 

The sample includes several very large and luminous 
members with well-extended rotation curves, e.g., UGC 
2885, NGC 801, and NGC 2903. They have high sur- 
face brightness (HSB), massive stellar component, and 
low gas content. Typically, their rotation curve rises 
steeply to a maximum and declines slowly into an al- 
most horizontal asymptote. There are also a number of 
dwarf, gas-dominated, and low-surface brightness (LSB) 
galaxies, e.g. DDO 168. There is no conspicuous max- 
imum, and in some galaxies not even a flat asymptote 
to their rotation curve. It is generally believed that de- 
viations from the classical dynamics is more pronounced 



in LSBs than in HSBs. (McGaugh fc de Blok 1998a 



ISanders fc Noordermeer 20071 IGentile et al. 20Tol) 

Twenty-eight members of the sample, of both HSB 
and LSB types, arc located in the Ursa Major cluster of 
galaxies, believed to be at the distance of about 15.5 Mpc 
d Tally fc Verheijen 1997 1. Seven of the galaxies, listed in 
Table HI have central bulges and are treated differently, 
the reason is explained below. For a full description of the 
sample the interested reader is referred to Sanders and 
McGaugh (2002). 



4. Constructing rotation curves 

We calculate the rotation speed of a test object circling the 
galaxy as a function of distance from the galactic center 
and the distribution of the detectable matter in the galaxy. 
The procedure we follow is almost that of Sanders and 
McGaugh (2002): 

In order to calculate the MOND rotation curves, the 
first step is determining the Newtonian acceleration of the 
detectable matter, via the classical Poisson equation. 
Given the Newtonian acceleration, the effective accelera- 
tion is calculated from the MOND Eq. For MOG, we 
approximate the galaxy by a spherically symmetric sys- 
tem. The error committed in this approximation, as de- 
scribed in Binney & Tremaine (1987), is on the order of 
15%. 



We assume a constant M^/L ratio throughout the 
galaxy, though this is not strictly the case, because of the 
color gradient in spiral galaxies. However, in seven bulged 
spirals we find assigning different M^,/L to the bulge and 
the disk improves the fit to the observed data. 

We assume the HI gas is in co-planer rota- 
tion about the center of the galaxy, an assumption 
which may not hold in galaxies with strong bars 
dSanders fc McGaugh 20021 ). 

Given the observed distribution of the baryonic mat- 
ter (stellar and gaseous disks, plus a spheroidal bulge, if 
present), the effective radial gravitational force, and sub- 
sequently the circular speed, is calculated from Eqs. ([Ij 
and (O - (inj. Fitting of the calculated rotation curves 
to the observed data points is achieved by adjusting the 
A'h/L ratio, through a least-square x^, defined as 



-E 

N ^ 



N 



theory ohs) 



N 



(9) 



i=l 



where ai is the observational uncertainty in the rotation 
speeds. The M^,/ L ratio of the disk and of the bulge are 
our ultimate results. 



5. Numerical results 

All models trace the observed data with reasonable de- 
tail. The best-fit and M,/L values are listed in Table 
[TJ Figures [T]- El show fits of theoretically constructed rota- 
tion curves to the observations of 46 galaxies. The general 
trend of HSB curves ( steep rise to a maximum followed 
by gradual decline to an almost flat asymptote), and of 
LSB curves ( slow rise often with no asymptote in sight) 
are evident. 

Seven galaxies have prominent bulge components. One 
expects a bulge with an older population of stars to have a 
higher Af*/L ratio than a disk with a younger population. 
Therefore, to obtain a better fit for these galaxies, we have 
allowed the model to choose different M*/Ls for the disk 
and the bulge. The result, shown in Tabled confirms the 
expectation. The minimum y^s of Table [21 are detectably 
lower than those of the corresponding entries of Table [H 
obtained by a single M^jL fit. Nevertheless, one galaxy in 
M0ND3, NGC 801, and four in MONDl, NGC 801, NGC 
5371, UGC 2885, and NGC 5907, predict untenably lower 
M^jL for the bulge than for the disk. 

McGaugh (private communication) advises us that this 
oddity might partially be owing to the low resolution of 
HI data and/or the sharp rise of v(r) in the bulge. If v(t) 
is not quite resolved, one tends to underestimate it, and 
the M^jL of the bulge with it. In many of these galax- 
ies it is indeed not obvious if the inner component is re- 
ally a "bulge" in the classical sense of a 3D component 
with an r^l^ profile. It is also said that neglecting the 
flattening of a bulge leads to an error in its mass-to-light 
ratio (Noordermeer 2008). In NGC 6946, there is a tiny 
bulge (4% of the total light in B-band). It shows up as a 
small kinematic bump in the inner 1 kpc in high-resolution 
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Fabry-Perot data (jBlais-Ouellette et al. 2004p . Moreover, 
the HI distribution in NGC 6946 is not symmetric in the 
galactic plane. It is patchy and seems to deviate from cir- 



cular orbits ( Carginan et al. 1990 1 



Ten entries in Table [T] have unacceptably large x^s. 
Suspecting the failure of the assumption of constant 
M^.IL^ we have followed Barnes et al. (2007), and ex- 
amined a radially varying, M^jL = [M^,/ L)q -\- mr. The 
best-fit values of the constants {M^,/L)o and m for some 
of them are displayed in Table |31 The slope, to, is much 
too small to result in appreciably lower x^- The reason for 
the failure should lie elsewhere. For example, the assump- 
tion of cold unobservable molecular gas in the galactic disk 
(Tiret & Combes 2009) leads to better fits with lower x^- 



6. Color - A'U/L correlation 

How realistic are the inferred M^,/L ratios? Stellar pop- 
ulation synthesis (SPS) models predict a linear rela- 
tion between colors and M^/L ratios. Redder galaxies 

[Bell fc de Jong 2001[ 



should have higher M^/L ( see, e.g.. 



IBell et al. 2"003l IPortinari et al. 2004p . The slope of this 
linear relation does not depend on the exact details of the 
history of star formation, i.e. the assumed IMF. But de- 
pending on how many stars are present at the low-mass 
end of the stellar IMF, the color-M*/L curve shifts up 
and down. This is because low-mass stars contribute sig- 
nificantly to the mass of a population, but not as much to 



its luminosity and color ( Bell fc de Jong 2001 ). 

In the SPS scheme Salpeter's (1955) IMF overesti- 
mates the Mif/L ratios of many of the galaxies and vi- 
olates the condition of 'less disk mass than the mass of 
maximum disk'. To remedy this. Bell et al. (2003) scale 
down Salpeter's IMF and come up with a limit for the 
color- A/* /L relation above which the physical viability is 
not guaranteed. Their suggested relation is 



log(M,/Ls) = 1.74(5 -V)- 0.94. 



(10) 



There are other IMF that lead to slightly different re- 
lations. For example, based on an analysis of the verti- 
cal velocity dispersion of stars, Bottema (1997) argues 
for a substantially submaximal M/ L ratio for all disk- 
dominated galaxies. Alternatively Kroupa (2001) intro- 
duces a turnover at the low-mass end of his IMF. 

In Eq. (ITUl) the slope 1.74 is not sensitive to variations 
in IMF, but the y-intercept is. To obtain the equivalent 
relation for standard Salpeter's, Kroupa's, and Bottema's 
IMF one should shift Eq. and the plots in Fig. (g]) up 
and down by about (0.15, -0.15, -0.35) dex, respectively 
(jBell et al. 20031) . 

In Fig. m we contrast M^jL ratios of the three gravity 
models against the predictions of SPS, where we use the 
B-band luminosities of Sanders & McGaugh (2002). In 
each frame the solid line is the best fit to the data points 
obtained from the analysis of the rotation curves. The 
theoretical SPS predictions of Bell & de Jong (2001), and 
Bell et al. (2003) for different IMFs are also plotted. The 



slope of MONDl, 1.78 ± 0.23, of M0ND2, 1.81 ± 0.21, 
and of M0ND3, 1.75 ± 0.26 are reasonably close to that 
of Eq. pUj) . The corresponding j/-intercepts, —0.88 ±0.14, 
— 1.06 ± 0.12, and —1.13 ± 0.15, respectively, are also in 
harmony with that of Eq. ([T0|) . The uncertainties in slopes 
and y-intercepts are in the 1-cr error range. The errors in 
^/-intercepts are small enough to enable one to distinguish 
between different IMFs. MONDl falls somewhere between 
standard Salpeter's and scaled Salpeter's IMF. M0ND2 
agrees with Kroupa's and scaled Salpeter's IMF. M0ND3 
is in good agreement with Kroupa's IMF. 

The slope for MOG, 1.06 ± 0.21, cannot be reconciled 
with SPS predictions. The case of NEW is also question- 
able. Although the slope, 2.33 ± 0.67, is consistent with 
1.74, considering its large error bar, the dispersion of the 
simulated data points is too large to conclude a meaning- 
ful color- M*/i correlation. 

Any alternative gravity can have a dark matter equiv- 
alent. Deviations from the Newtonian gravity can be at- 
tributed to a hypothetical dark entity and a dark den- 
sity profile calculated through Poisson's equation, for in- 
stance. One feature however distinguishes this interpreta- 
tion from the conventional CDM scenarios. Here, there is 
a well defined relation between the baryonic matter and 
its so-interpreted dark companion. While in CDM models, 
baryonic and dark matters may co-exist independently. In 
our opinion, the reason for the good agreement of MOND 
with SPS predictions and non-compliance of CDM with 
it lies in the existence or non-existence of this relation 
between the observable and non-observable matters. In 
MOND, baryonic matter plays a pivotal role, and the dark 
entity owes its existence to it. This is not the case in CDM. 
Dark matter is allowed to play a role independently ofthe 
observable matter. As for MOG, we are not in a position 
to express an opinion. 

Let us summarize our conclusion: a) the SPS scheme 
can distinguish between different gravity models and b) 
the two together can choose between different IMFs. The 
mere fact that a gravity theory reproduces the observed 
rotation curves satisfactorily does not tell the whole story. 



7. Concluding remarks 

At least on galactic scales, dynamics of spirals cast doubt 
on the viability of the classical theories of gravitation. A 
number of alternative theories are capable of reproducing 
the rotation curves of spirals with acceptable detail, a non- 
trivial fact that deserves attention. In this paper we used 
two alternative theories of gravitation, MOND and MOG, 
and a CDM model to deduce the dynamics of a well-sized 
sample of high- and low-surface brightness galaxy types, 
and checked the results against observations. The models 
are not equivalent although they all simulate the rotation 
curves in more or less to the same degree of accuracy, . 

In MOND and MOG, rotation curves are constructed 
with only one free adjustable parameter, the stellar mass- 
to-light ratio. This is in contrast to the CDM model, where 
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an additional parameter is needed to describe the dark 
component. 

There are cases of bulged galaxies where fits to obser- 
vations lead to lower A/*/L ratios for the bulge than for 
the disk. This might be owing to the low resolution of the 
HI data and to the inner v{r), and/or uncertain size of the 
bulge. 

Stellar population synthesis models impose constraints 
on M^,/L. Redder galaxies should have higher Af*/L ra- 
tios. MOND with different interpolating functions meets 
this expectation, albeit with different IMFs. This is re- 
markable, because there is no explicit /implicit connection 
between the basic tenets of the SPS and MOND. On the 
other hand, MOG does not meet the SPS constraints, and 
the data points of CDM show a large scatter, preventing 
one to conclude a meaningful color- M*/L correlation. 

The SPS predictions of M*/L ratios are sensitive 
to the adopted IMF. The M*/L ratios inferred from 
Salpeter's IMF are notably higher than those obtained 
from Kroupa's. M0ND3 favors Kroupa's IMF. It produces 
lower M*/_L and implies lesser disk masses. 
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Galaxy (Type) B-V JrJI ^ JrJI ^ ITJL ^ ATjI ^ ATjI 

Mondl Mondl Mond2 Mond2 MondS Mond3 MOG MOG NFW NFW 



HSB Galaxies 



M 33 (Sc) 


0.55 


0.6 


30.22 


0.4 


36.47 


0.3 


39.25 


0.8 


45.03 


0.2 


125.00 


NGC 300 (Sc) 


0.58 


0.7 


2.26 


0.5 


2.51 


0.4 


2.35 


1.3 


2.54 


0.1 


2.94 


NGC 2903 (Sc) 


0.55 


3.0 


6.07 


2.2 


5.72 


1.7 


4.79 


2.4 


7.18 


2.4 


3.60 


NGC 3726 (SBc) 


0.45 


1.0 


3.47 


0.7 


3.86 


0.5 


4.60 


0.9 


5.00 


0.6 


3.03 


NGC 3769 (SBb) 


0.64 


1.2 


0.75 


0.9 


0.65 


0.7 


0.70 


1.4 


1.06 


0.7 


0.42 


NGC 3877 (Sc) 


0.68 


1.7 


2.66 


1.2 


2.62 


0.9 


2.60 


1.5 


3.07 


1.1 


2.51 


NGC 3893 (Sc) 


0.56 


1.7 


4.12 


1.3 


2.11 


1.0 


1.88 


1.6 


2.76 


1.2 


1.89 


NGC 3949 (Sbc) 


0.39 


0.8 


5.34 


0.6 


3.88 


0.5 


3.97 


0.8 


5.00 


0.5 


3.29 


NGC 3953 (SBbc) 


0.71 


2.7 


1.13 


2.0 


0.48 


1.5 


0.47 


2.2 


0.44 


2.1 


1.05 


NGC 3972 (Sbc) 


0.55 


1.5 


3.22 


1.0 


3.03 


0.8 


3.08 


1.6 


2.85 


0.2 


1.25 


NGC 3992 (SBbc) 


0.72 


4.9 


0.65 


3.6 


0.88 


2.7 


1.02 


3.6 


2.08 


4.5 


4.11 


NGC 4013 (sb) 


0.83 


3.1 


1.37 


2.3 


1.62 


1.8 


2.24 


2.7 


2.05 


2.2 


1.63 


NGC 4051 (SBbc) 


0.62 


1.2 


0.88 


0.9 


0.78 


0.7 


0.90 


1.1 


0.76 


0.7 


0.78 


NGC 4085 (Sc) 


0.47 


1.1 


6.84 


0.8 


6.02 


0.6 


6.03 


1.1 


6.94 


0.5 


3.93 


NGC 4088 (Sbc) 


0.51 


1.1 


1.49 


0.8 


1.62 


0.6 


1.75 


1.0 


1.97 


0.8 


1.24 


NGC 4100 (Sbc) 


0.63 


2.4 


2.07 


1.7 


1.76 


1.3 


1.58 


2.0 


1.84 


1.9 


1.69 


NGC 4138 (Sa) 


- 


3.5 


1.61 


2.7 


0.98 


2.0 


0.81 


3.2 


1.20 


2.9 


1.09 


NGC 4157 (Sb) 


0.66 


2.4 


0.92 


1.7 


0.85 


1.3 


0.87 


2.0 


0.84 


1.7 


0.89 


NGC 4217 (Sb) 


0.77 


2.2 


3.95 


1.6 


2.92 


1.2 


2.90 


1.9 


2.63 


1.5 


3.02 


NGC 4389 (SBbc) 


- 


0.4 


5.36 


0.3 


5.42 


0.2 


5.57 


0.6 


6.33 


0.1 


4.33 


NGC 5585 (SBcd) 


0.46 


0.5 


10.43 


0.4 


10.22 


0.3 


10.46 


1.1 


16.87 


0.1 


19.03 


NGC 6946 (SABcd) 


0.40 


0.5 


11.46 


0.4 


17.90 


0.3 


21.61 


0.5 


31.75 


0.3 


1.01 


NGC 7793 (Scd) 


0.63 


1.2 


1.48 


0.9 


1.44 


0.6 


1.48 


1.5 


1.02 


0.8 


1.91 


UGC 6399 (Sm) 




1.0 


0.16 


0.8 


0.17 


0.6 


0.10 


1.8 


0.04 


0.1 


1.48 


UGC 6973 (Sab) 
NGC 801' (Sc) 
NGC 2998' (SBc) 




2.7 


20.46 


2.2 


11.81 


1.7 


10.57 


2.6 


20.24 


2.0 


6.48 


0.61 


1.2 


23.14 


1.0 


14.75 


0.8 


14.15 


1.2 


23.90 


1.2 


49.61 


0.45 


1.2 


2.64 


0.9 


2.43 


0.7 


2.96 


1.0 


2.35 


1.4 


6.80 


NGC 5371' (S(B)b) 


0.65 


1.6 


10.02 


1.2 


8.32 


0.9 


6.93 


1.3 


6.65 


1.5 


7.65 


NGC 5533' (Sab) 

NGC 5907' (Sc) 


0.77 


3.3 


2.33 


2.6 


1.61 


2.1 


1.11 


3.8 


8.50 


5.6 


19.61 


0.78 


4.0 


2.93 


2.8 


3.82 


2.1 


4.27 


3.0 


6.10 


3.5 


11.23 


NGC 6674' (SBb) 


0.57 


2.7 


10.96 


2.0 


7.65 


1.6 


6.64 


2.6 


41.06 


4.1 


66.96 


UGC 2885' (Sbc) 


0.47 


1.5 


2.80 


1.2 


2.98 


0.9 


3.04 


1.4 


6.64 


1.9 


14.85 


LSB Galaxies 


DDO 168 (SO) 


0.32 


0.2 


11.50 


0.1 


14.35 


0.1 


21.56 


1.5 


14.64 


0.1 


26.67 


NGC 247 (SBc) 


0.54 


1.1 


3.71 


0.8 


3.91 


0.7 


4.16 


2.0 


3.74 


0.1 


10.34 


NGC 1560 (Sd) 


0.57 


1.1 


3.35 


0.6 


1.94 


0.3 


1.52 


4.6 


10.56 


0.1 


17.00 


NGC 3917 (Scd) 


0.60 


1.3 


4.49 


0.9 


4.58 


0.7 


4.58 


1.4 


4.03 


0.2 


6.51 


NGC 4010 (SBd) 


0.54 


1.4 


1.81 


1.0 


1.74 


0.8 


1.76 


1.7 


1.16 


0.1 


2.42 


NGC 4183 (Sa) 


0.39 


0.7 


0.98 


0.5 


0.98 


0.4 


1.09 


1.0 


1.54 


0.4 


0.20 


UGC 128 (Sdm) 


0.60 


1.1 


0.48 


0.8 


0.54 


0.6 


0.63 


1.9 


0.36 


0.1 


2.57 


UGC 6446 (Sd) 


0.39 


0.5 


2.30 


0.4 


3.29 


0.3 


4.49 


1.2 


2.35 


0.1 


0.30 


UGC 6667 (Scd) 


0.65 


1.0 


0.94 


0.8 


0.88 


0.6 


0.69 


1.9 


0.59 


0.1 


3.95 


UGC 6917 (SBd) 


0.53 


1.4 


0.64 


1.0 


0.69 


0.8 


0.72 


2.0 


0.84 


0.1 


0.49 


UGC 6923 (Sdm) 




0.8 


1.17 


0.6 


1.16 


0.4 


1.03 


1.4 


2.28 


0.1 


0.56 


UGC 6930 (SBd) 


0.59 


0.8 


0.28 


0.6 


0.34 


0.4 


0.54 


1.2 


0.34 


0.2 


0.19 


UGC 6983 (SBcd) 


0.45 


1.7 


1.30 


1.2 


1.46 


0.9 


1.68 


2.3 


1.90 


1.1 


0.54 


UGC 7089 (Sdm) 




0.2 


0.14 


0.2 


0.40 


0.1 


0.25 


0.6 


0.11 


0.9 


3.18 



Table 1. Best-fit reduced and M^/L values of 32 HSB and 14 LSB galaxies in MOND with different interpo- 
lating functions, MOG. and NFW models. For explanation of models see text. Hubble types are from NASA/IPAC 
Extragalactic Database (NED). Bulged galaxies are marked by a superscript 
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Galaxy 


XMOND3 


/ M* \disk 
\~ )MOND3 


/ M* \bulge 

\ L )mOND3 


Xmondi 


/ M* \disk 

J MONDl 


( M* \bulge 

\ L IMONDI 


NGC 801 


12.69 


i.r 


0.7 


17.72 


2.2* 


1.1 


NGC 2998 


2.22 


0.6 


0.8 


2.50 


1.2 


1.3 


NGC 5371 


6.34 


0.8 


1.0 


9.96 


1.7* 


1.6 


NGC 5533 


0.98 


0.7 


2.2 


2.12 


0.1 


3.7 


NGC 5907 


3.5 


1.8 


3.2 


2.80 


4.1* 


3.6 


NGC 6674 


6.32 


0.4 


1.8 


10.83 


1.5 


2.9 


UGC 2885 


2.07 


0.8 


1.1 


2.80 


1.5* 


1.4 



Table 2. Same as Table 1 for seven bulged galaxies. Disks and bulges have different M*/L ratios. One galaxy in 
M0ND3 model and four in MONDl's predict untenably lower M*/L for the bulge than for the disk. These are marked 
by asterisks in the M*/L columns. 

Galaxy Xmond3 {^)o,mond3 niMONDS Xaiondi {^)o,mondi mMONPi 

M 33 37.4 0.4 -0.01 27 0.6 0.01 

DDO 168 20.19 0.1 -0.02 10 0.1 0.07 

NGC 801 14.15 0.8 0.00 19.88 1.1 0.01 



Table 3. Linear, M*/L = (M*/L)o + mr, fit for 3 cases not well explained by either M0ND3 formalism or by that 
of MONDl. A comparison with constant M*/L of these galaxies (Table 1) shows no significant improvement. 
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Fig. 1. Rotation curves of 32 mainly HSB galaxies. The points with vertical error bars are the observed data. The 
dotted (black) and short dashed lines are contributions of the gaseous and stellar components to the Newtonian rotation 
speeds, respectively. The long dashed (blue) line is that of MONDl. The solid line is the rotation curve constructed 
through M0ND3. The M*/Ls of M0ND3 are used in plotting the stellar component. The dashed-dotted (green) and 
dotted (red) lines are those of the MOG and NFW models, respectively. 
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Fig. 2. Fig. 1 continued. The last seven galaxies have a bulge component, depicted as dashed-dotted lines. 
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DDO 1 68 NGC 247 





Fig. 3. Rotation curves of 14 LSB galaxies. Legend as in Fig. 1. 
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Kroupa IMF 
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Scaled Salpeter IMF 
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MOG 





Kroupa IMF 
Salpeter IMF 
1.06(B-V)-0.37 
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Fig. 4. Plots of Mif/L versus B — V . The solid line in each panel is the best fit to the prediction of the model in ques- 
tion. The slopes and y-intercepts of best-fitted lines are shown in the panels. Slopes of MONDl, M0ND2, and M0ND3 
are reasonably close to the prediction of SPS, 1.74. That of MOG is not. Data points in NFW panel are much dispersed 
to conclude a meaningful correlation. In MONDl, M0ND2, and M0ND3 errors in the slopes, ±(0.23,0.21,0.26), re- 
spectively, are small enough to distinguish one model from the other. The errors in the j/-intercepts, ±(0.12, 0.14, 0.15), 
respectively, are small enough to distinguish one IMF from the other. The other lines are the theoretical predictions 
of SPS with different IMFs ( |Bell fc de Jong 200T|IBell et al. 20 031). They are included here for comparison. They have 
almost the same slope, but different j/-intercepts. 



